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Performing a shell model calculation for heavy nuclei has been a long-standing problem in nuclear 
physics. Here we propose one possible solution. The central idea of this proposal is to take the 
advantages of two existing models, the Projected Shell Model (PSM) and the Fermion Dynamical 
Symmetry Model (FDSM), to construct a multi-shell shell model. The PSM is an efficient method of 
coupling quasi-particle excitations to the high-spin rotational motion, whereas the FDSM contains 
a successful truncation scheme for the low-spin collective modes from the spherical to the well- 
deformed region. The new shell model is expected to describe simultaneously the single-particle and 
the low-lying collective excitations of all known types, yet keeping the model space tractable even 
for the heaviest nuclear systems. 

PACS numbers: 21.60.Cs 



I. INTRODUCTION 

Except for a few nuclei lying in the vicinity of shell 
closures, most of the heavy nuclei are difficult to de- 
scribe in a spherical shell model framework because of 
the unavoidable problem of dimension explosion. There- 
fore, the study of nuclear structure in heavy nuclei has 
relied mainly on the mean- held approximations, in which 
the concept of spontaneous symmetry breaking is applied 
0,0- However, there has been an increasing number of 
compelling evidences indicating that the nuclear many- 
body correlations are important. Thus, the necessity of a 
proper quantum mechanical treatment for nuclear states 
has been growing, and we are facing the challenge of un- 
derstanding the nuclear structure by going beyond the 
mean-held approximations. 

Demand for a proper shell model treatment arises also 
from the nuclear astrophysics. Since heavy elements are 
made in stellar evolution and explosions, nuclear physics, 
and in particular nuclear structure far from stability, en- 
ters into the stellar modeling in a crucial way. The nu- 
cleosynthesis and the correlated energy generation is not 
completely understood, and the origin of elements in the 
cosmos remains one of the biggest unsolved physics puz- 
zles. Nuclear shell models can generate well-defined wave 
functions in the laboratory frame, allowing us to com- 
pute, without further approximations as often assumed 
in the mean-field approaches, the quantities such as tran- 
sition probabilities, spectroscopic factors, and /?-decay 
rates. These quantities provide valuable structure in- 
formation to nuclear astrophysics. In fact, the nuclear 
shell model calculations could strongly modify the results 
of nuclear astrophysics, as the recent work of Langanke 
and Martinez-Pinedo has demonstrated (see, for exam- 
ple, Ref. 0). 

Tremendous efforts have been devoted to extending the 
shell model capacity from its traditional territory of the 
SG?-shell to heavier shells. Over the years, one has looked 
for possible solutions in the following two major direc- 



tions. In the first direction, one employs rapidly growing 
computer power and sophisticated diagonalization algo- 
rithms to improve the traditional shell model code. The 
shell model code ANTOINE 4j is a representative exam- 
ple of the recent developments along this line. Using this 
code the deformed N Z nuclei up to mass A ~ 50 can 
be well explained. The recent record example performed 
by this code is the full fp shell calculation of A = 52 
nuclei [5j, with the basis dimensions in excess of 100 mil- 
lion. 

While in principle, it does not matter how to prepare 
a shell model basis, it is crucial in practice to use the 
most efficient one. Moreover, feasibility in computation 
is not our only concern. The other important aspect 
of using an efficient basis is that it may have a good 
classification scheme, such that a simple configuration in 
that basis corresponds approximately to a real mode of 
excitation. This not only can simplify the calculations, 
but also make the physical interpretations of results more 
easy and transparent. 

In the second direction, which is defined in a much 
wider scope, one employs various methods in seeking ju- 
dicious truncation schemes. Such schemes should con- 
tain the most significant configurations, each of which 
can be a complicated combination in terms of the orig- 
inal shell model basis states. In this way, the basis di- 
mension can be significantly reduced and the final diag- 
onalization is performed in a much smaller space, thus 
making a shell model calculation for heavy nucleipossi- 
ble. The early MONSTER- VAMPIRE approach and 
the recent Monte Carlo Shell Model 0, |g| are examples 
along this line. Nevertheless, numerical calculations re- 
quired by these models are still quite heavy, which may 
make a systematical application difficult. 

There are two other existing models that belongto the 
second category: the Project Shell Model (PSM) Q, and 
the Fermion Dynamical Symmetry Model (FDSM) [hJ. 
In the PSM, the shell model basis is constructed by choos- 
ing a few quasi-particle (qp) orbitals near the Fermi sur- 
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faces and performing angular-momentum and particle- 
number projection on the chosen configurations. By tak- 
ing multi-qp states as the building blocks, the PSM has 
been designed to describe the rotational bands built upon 
qp excitations |9(. The PSM has been rather successful 
in calculating the high-spin states of well-deformed and 
superdeformed nuclei. For lighter nuclei where the large- 
scale shell model calculation is feasible studies for 
the deformed 48 Cr [l^ and the superdeformed 36 Ar ^| 
have demonstrated that the PSM calculation can achieve 
a similar accuracy in describing the data. In the FDSM, 
the truncated basis is built by the symmetry detected S- 
and D-pairs, assuming that these are the relevant degrees 
of freedom for the low- lying collective motions. Having 
these pairs as the building blocks, the FDSM can pro- 
vide a unified description for the low-spin collective ex- 
citations from the spherical to the well-deformed region 



It is clear that the PSM and the FDSM follow the shell 
model philosophy and both have their own shell model 
truncation scheme. However, the truncations emphasize 
different excitation modes, which are contained in one 
model but are absent in the other. An idea emerges 
naturally that one may combine the advantages of the 
two models to construct a new shell model for heavy nu- 
clei. The question is how. The PSM is a microscopic 
approach employing the deformed intrinsic states and 
the projection method, while the FDSM is a fermionic 
model based on the group theory. The crucial step that 
leads us to connect these two different approaches is made 
through the recent recognition |14i that the numerical 
results obtained by the PSM exhibit, up to high angular 
momenta and excitations, a remarkable one-to-one cor- 
respondence with the analytical SU(3) spectrum of the 
FDSM. This suggests that the projected deformed-BCS 
vacuum has at the microscopic level SU(3)-\ike struc- 
tures which are very close to the representations of the 
SU(3) dynamical symmetry of an S-D fcrmion-pair sys- 
tem. This recognition has motivated us to propose a 
multi-shell shell model for heavy nuclei. Hereafter, we 
shall call it Heavy Shell Model, or HSM for short. 



In the following section, the PSM and the FDSM will 
be briefly reviewed. The emphasis will be given on 
the discussion of the advantages and deficiencies of each 
model. In Section 3, the connection between the two 
models will be explored. In Section 4, we will discuss 
in detail how the two models are integrated to form the 
HSM. We will give the basis states and the basis trunca- 
tions for the well-deformed, transitional, and spherical re- 
gions. The effective interactions and the general method 
for evaluating the projected matrix elements will also be 
discussed in this section. Finally, the paper will be sum- 
marized in Section 5. 



II. PROJECTED SHELL MODEL AND 
FERMION DYNAMICAL SYMMETRY MODEL 

In this section, we introduce the basic structure of the 
PSM and the FDSM. For each model, we point out the 
main features and the limitations. For interested readers, 
we refer to the review article of the PSM and of the 
FDSM [nj. 



A. Projected Shell Model 

The PSM begins with the deformed (e.g. the Nilsson- 
type) single particle basis, with pairing correlations in- 
corporated into the basis by a BCS calculation for the 
Nilsson states. The basis truncation is first implemented 
in the multi-qp states with respect to the deformed- 
BCS vacuum |$); then the angular- momentum and the 
particle-number projection are performed on the selected 
qp basis to form a shell model space in the laboratory 
frame; finally a shell model Hamiltonian is diagonalized 
in this projected space. 

If aj, and a\ are the qp creation operators, with the in- 
dex Ui (jri) denoting the neutron (proton) quantum num- 
bers and running over properly selected single-qp states, 
the multi-qp bases of the PSM are given as follows: 

even-even nucleus: {|$) , 4*4- |$) , 4 4 4j |$) , 

odd-i/ nucleus: {a+. |$> , 4,4,4* 1$) , ' ' • K 1 ) 
odd-7r nucleus: {4 ; |$) > a t i a l a l k 1^) >" " } 
odd-odd nucleus: {a) u .a\. |$) , ■ • ■ } 

In bases (1), "• • • " denotes those configurations that con- 
tain more than two like-nucleon quasi-particles. If one is 
interested in the low-lying states only, they can practi- 
cally be ignored because these configurations have higher 
excitation energies due to mutual blocking of levels. The 
bases (1) can be easily enlarged by including higher or- 
der of multi-qp states if necessary. If the configurations 
denoted by "• • • " are completely included, one recovers 
the full shell model space written in the representation 
of qp excitation. 

In the qp basis, the truncation can now be easily im- 
plemented by simply excluding the states with higher 
energies. Usually, only a few orbitals around the Fermi 
surfaces are sufficient for a description of the low- lying qp 
excitations. The truncation is thus so efficient that di- 
mension never poses a problem even for superdeformed, 
heavy nuclei. 

After the truncation is implemented in the multi-qp 
basis, the shell model space can be constructed by the 
projection technique 0: 

\qKIM) = P{f K |$ 9 ) with P I M N K = P I MK P N : (2) 
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where denotes the qp-basis given in (1) with q mean- 
ing the multi-qp configuration, and 

p mk = J dnDl IK (n)R(n) 

= J_ f dfc-Ufi-*)* ( 3 ) 
2tt J 

are the angular momentum and particle number pro- 
jection operators, respectively. In Eq. ©, D I MK is the 
-D-function [l^, R the rotation operator, fl the solid 
angle, N the number operator, and <f> the gauge an- 
gle. If one keeps the axial symmetry in the deformed 
basis, D\ 1K in Eq. J3J reduces to the small d-function 
and the three dimensions in Q reduce to one. The 
eigen-values E and the corresponding wave functions 

|*m) = J2 q K F qK\<l KIM ) are tnen obtained by solv- 
ing the following eigen-value equation: 



£ i H U </*> - E N U } K'k> = o, (4) 

qK 

where H qK q , K , and N* K q , K , are respectively the matrix 
elements of the Hamiltonian and the norm 

HU q >K> = (qKI\H\q'K'I) = (4 q \H P^,|$ q ,) 

(5) 

N} Kq , KI = (qKI\q'K'I) = (<Pq\Pk N K >\<S> q >)- 

The PSM uses a large size of single-particle (s.p.) 
space, which ensures that the collective motion is de- 
fined microscopically by accommodating a sufficiently 
large number of active nucleons. It usually includes three 
(four) major shells each for neutrons and protons in a cal- 
culation for deformed (superdeformed) nuclei. The effec- 
tive interactions employed in the PSM are the separable 
forces. The Hamiltonian takes the following form: 

. . , (6) 

The first term Hq in of Eq. J^J is the spheri- 
cal single-particle Hamiltonian and the remaining terms 
are residual quadrupole-quadrupole, monopole-pairing, 
and quadrupole-pairing interactions, respectively. The 
strength of the quadrupole-quadrupole force is deter- 
mined in a self-consistent way that it would give the 
empirical deformation as predicted in the variation cal- 
culations. The monopole-pairing strength is taken as the 
form Gm = G j A (A is the mass number) with G being 
adjusted to yield the known odd-even mass differences. 
The quadrupole-pairing strength Gq is assumed to be 
about 20% of G M @. 

The one-body operators (for each kind of nucleons) in 



Eq. JBJ) are of the standard form 

a, a' 

P' = ^£44 (7) 

a 
a, a' 

where Q MQQ < = (a\Q2fi\cx.') is the one-body matrix ele- 
ment of the quadrupole operator and c' a the nucleon cre- 
ation operator, with a standing for the quantum numbers 
of a s.p. state in the spherical basis (a = {n£jm}). The 
time reversal of c\ is defined as c a = (— ) J_m c„^'_ m . 

For a description of rotational bands associated with 
a well-deformed minimum, the PSM is a highly efficient 
truncation scheme. Diagonalization for a heavy nucleus 
can be done almost instantly, yet the results are often sat- 
isfactory. The reason for the success is because the major 
part of pairing and quadrupole correlations has already 
been built in the basis through the use of deformed basis 
and the BCS formalism. Therefore, a small configuration 
space with a few s. p. orbits around the Fermi surface can 
already span a very good basis for the low-lying excita- 
tions. Note that each of the configurations in the basis 
is a complex mixture of multi-shell configurations of the 
spherical shell model space. Although the final dimen- 
sion of the PSM is small, it is huge in terms of original 
shell model configurations. In this sense, the PSM is a 
shell model in a truncated multi-major-shcll space. 

The pleasant features of the PSM make it a frequently- 
used model in the high-spin physics. Many applications 
can be found in the review article 9]. The recent publi- 
cations include the study of superdeformed structure in 
a wide range of different mass regions Hjl [lSl fl9l l20| , 
the study of the origin of identical bands [21j and of the 
high-K states [H. 

While the PSM is an efficient shell model for deformed 
systems with rotational behavior, it becomes less valid 
when going to the transitional region, and eventually 
loses its applicability for spherical nuclei. Moreover, it 
can not efficiently describe the j3- and 7-vibrations. Al- 
though such collective modes can in principle be obtained 
by mixing a large amount of excited qp-configurations, it 
suffers in practice from a similar dimension problem as 
in conventional shell models. The main reason for these 
shortcomings is due to the use of the simple BCS- vacuum 
which contains only the properties of the ground-state ro- 
tational band but not those of the collective vibrations. 



B. Fermion Dynamical Symmetry Model 

If we say that the PSM is a shell model in a truncated 
multi-major-shell space, then the FDSM is a shell model 
in a truncated one-major-shell space. The truncation is 
based on the consideration that the like-nucleons prefer 
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to form coherent 5* (angular-momentum L = 0) and D 
(L = 2) pairs. One may thus assume that a closed sub- 
space built up by these S-D pairs is mainly responsible 
for nuclear low-lying collective motions. Such an S-D 
subspace can be carved out by a symmetry requirement 
that the S and D creation (annihilation) operators to- 
gether with necessary minimum amount of number con- 
serving operators form a closed Lie algebra. To this end, 
a k-i basis is introduced (see Ref. [T^j an d the references 
cited therein) 



E 



-t 



(8) 



where k (pseudo-orbital angular-momentum) and i 
(pseudo-spin) could be either k = 1 and i — any half 
integer, or % — 3/2 and k — any integer. This basis must 
uniquely reproduce the normal-parity levels |23| in that 
shell by k-i angular momentum coupling, no more and 
no less. With this k-i basis, the coherent S and D pairs 
and the multipole operators P rf j,, which are necessary in 
order to form a closed Lie algebra, are found to be as 
follows H: 



(for any k and i) 
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(fc-active: k = 1, r < 2) (i-active: i = 3/2, r < 3) 

(11) 

where the symbol [ ] denotes angular-momentum 
coupling, and the time reversal is defined as 
bfcmkimj = {-) k ~ mk+t ~ mi b k - mk i- mi - The operator 
set {S\ S, DL D^, P rfJ ,} forms a closed Lie algebra of 
either Sp(6) (fc-active) or 5*0(8) (i-active), depending 
on the level structure of the valence shell 

Once the S-D subspace is carved out, the form of the 
effective Hamiltonian (restricted to a two-body interac- 
tion) in this truncated space is uniquely determined: 

H = U ° + H ^ 1 H ™ = ~E B 7 P r- P r 

(T — U.-K r 

H a = H£-G%S crJi S a -GZD aj -D< r -^2B?P?-P?(12) 

r>0 

It can be shown that the multipole operators P r with 
r = 0, 1 are proportional to the number operator in nor- 
mal parity levels n\ and the total angular-momentum 
J, respectively, while Pi is proportional to the effective 
quadrupole operator in the truncated space. The term 



Hq is a quadratic function of valence neutron and pro- 
ton numbers (ji*n v is included in the n-p interaction 
H U7r ). This Hamiltonian (Eq. |JT2J|) l° oks formally simi- 
lar to that of the PSM (Eq. ©) if we assume that only 
P>2 7^ in the summation. However, one should bear in 
mind that the FDSM Hamiltonian is written in a trun- 
cated one-major shell space and the s. p. energy splitting 
within the shell is neglected. This is the price the FDSM 
has to pay in order to meet the symmetry requirement 
so that solutions can be obtained with the aid of group 
theory. 

With this Hamiltonian, it can be shown that there 
exist analytical solutions in various dynamical symme- 
try limits, each one corresponding to a collective mode 
known experimentally: the SU(3)-\hmt in the fc-active 
shell and the S'0(6)-limit in the z-active shell correspond 
respectively to rigid and 7-soft rotors for well deformed 
nuclei, while the 5'C/(2)-limit in the /c-active shell and the 
S'0(5)-limit in the i-active shell correspond to a vibrator 
of spherical nuclei. Since the FDSM contains all major 
collective modes, the general feature of different collective 
motions arises naturally as the number of valence nucle- 
ons varies; namely, nuclei behave as a spherical vibrator 
near the closed shell and become a well-deformed rotor 
around the mid-shell. If the strengths of the Hamiltonian 
are properly chosen and the s. p. splitting is taken into 
account as a perturbation, the FDSM can even quantita- 
tively reproduce the low-lying spectra, B(E2)'s, ground 
state masses, etc., in a unified manner from the spherical 
to well-deformed region. 

Here, it may be appropriate to emphasize one remark- 
able result of the FDSM: There exists an one-to-one cor- 
respondence between the SU(3) irreps and the (3- and 
7- vibrations in deformed nuclei. As one can see from Ta- 
ble I, the /3-7 vibrations are microscopically the collec- 
tive D-pair excitations (the D^ ±1 excitations are forbid- 
den by the time reversal symmetry). The anharmonic 
behavior of the /3-7 vibrations is due to the finite parti- 
cle number effect. In the large n\ limit, ignoring 1/ni, 
the FDSM reproduces exactly the Particle-Rotor Model 
results. This means that the FDSM has discovered the 
relevant fermion degrees of freedom of nuclear collective 
motions, which are the S and D pairs. The S-D subspace 
is so compact that it never suffers from the dimension ex- 
plosion even for the heaviest nuclear systems. 

While with the S-D subspace the FDSM is able to 
provide a microscopic view to the low-lying collective mo- 
tions, it has difficulties in describing the qp excitations 
and the high-spin physics for lack of the s. p. degrees of 
freedom. In principle, this can be resolved by allowing 
a few pairs to break. The problem is that once the s. p. 
degrees of freedom open up, the dimension increases very 
rapidly. In addition, inclusion of s. p. degrees of freedom 
results in adding many new terms to the effective Hamil- 
tonian. To pin down so many coupling strengths in the 
Hamiltonian is a very difficult task, if not impossible. 
Therefore, even the number of broken pairs is limited to 
just one for proton and one for neutron, the model could 
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TABLE I: The SU(3) irreps in the FDSM and the collective f3-j vibrations. * 
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* The St/ (3) intrinsic states for each of the f3-y vibrational modes are listed in the first column. The second and the third 
column are the phonon excitation energies and the associated quantum numbers according to the Particle-Rotor Model. The 
4th and 5th column are the SU(3) irreps (A,/i) and the corresponding Young tablets [hi,h2,li3], respectively. The excitation 
energies obtained from the FDSM are listed in the last column as E>2AC(\, n), where AC(A, /i) is the changes of the 
expectation values of the 517(3) Casimir operator with respect to the ground state irrep C(m, 0), £>2 is the n-p quadrupole 
interaction strength, and hui = (3/2)ni_E>2. 



at best be applied to the fp shell nuclei, and could not 
go any further. 



III. CONNECTION BETWEEN PSM AND 
FDSM 

Let us summarize our main claims in the last section. 
Both the PSM and the FDSM are truncated shell model, 
aiming at grasping the essential ingredients to describe 
the low-lying physics. However, the emphasis in each 
model is different, which is reflected in their different 
truncation schemes. The PSM emphasizes the high-spin 
description of rotational states built upon qp-excitations 
associated with a well-deformed minimum, but it is not 
an efficient method for the collective vibrations. In con- 
trast, the FDSM has a well-defined classification for all 
known types of collective vibration, ranging from the 
spherical, via transitional, to the well-deformed region, 
but it is lack of the necessary degrees of freedom of qp- 
excitations. Thus, the main advantages of the PSM and 
the FDSM are mutually complementary to each other. 

If we could take the advantages of the PSM and the 
FDSM and combine them into a single model, the defi- 
ciencies of each model will be eliminated. At first glance, 
it is not obvious how to bridge these two different ap- 
proaches. In this section, we show that a realization 



of the combination is possible. The assertion is made 
based on the recent recognition that an SU (3) symmetry 
can emerge from the projected deformed-BCS vacuum. 
Having this as the basis, the full idea of the microscopic 
classification for the collective vibrations discovered by 
the FDSM could be adopted by the PSM while the lat- 
ter keeps its original features in building the shell model 
configuration space through projection. 



A. Emergence of SU(3) Symmetry in PSM 

It is remarkable [bH Il5| that the states numerically 
obtained by the PSM exhibit a beautiful one-to-one cor- 
respondence with the analytical SU(3) spectrum of the 
FDSM. To show this, the original PSM was extended in 
such a way that instead of a single BCS vacuum, the an- 
gular momentum projection was performed for separate 
neutron and proton BCS vacuum with the same defor- 
mation, and the projected states were coupled through 
diagonalization of the Hamiltonian The extension 

is necessary; Otherwise the PSM would have no collec- 
tively excited bands to compare with the FDSM. This 
procedure gives the usual rotational ground-band corre- 
sponding to a strongly coupled BCS condensate of neu- 
trons and protons, but also led to a new set of excited 
bands arising from the vibrations of relative orientation 
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FIG. 1: Comparison between the PSM (symbols) and the 
FDSM 517(3) limits (curves) assuming n^ ff = 32 and n? ff = 
16. The SU (3) quantum numbers (A, [i)k are marked for each 
band with degeneracy of k = K ma „ K m «, — 2, • • ■ , or 1, and 
K max = min{\, fj,}. The comparison takes the well-deformed 
nucleus 168 Er as an example. 



of the neutron and proton cores (the so-called "scissors" 
mode). 

As one can see from Fig. 1, the spectrum obtained from 
the PSM is, up to high angular momenta and excitation 
energies, nearly identical to that from the FDSM SU (3) 
formulism with (n° ff = 32, 0) and (n° ff = 16, 0) irreps, re- 
spectively assigned to the neutron and proton BCS vac- 
uum. The classification of the spectrum follows exactly 
the FDSM SU(3) reduction rules 
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1 ft, k + 1, • • • , A + (1 — k + 1, if k 7^ 0. 
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Not only the spectrum but also the B(E2)'s exhibit this 
correspondence 0. Note that there are many types of 
5(7(3) models, and that different 5(7(3) models have dif- 
ferent permissible irreps and reduction rules due to the 
different physical input, and therefore can lead to differ- 
ent band structures. Here, we emphasize that this SU(3) 
symmetry shown in Fig. 1 is of the FDSM type. This 
is because among the existing fermionic ST/ (3) models 
in nuclear physics, only the FDSM SU(3) formulism can 
naturally provide the required irreps and reproduce the 
band structure of the PSM. Other SU(3) models such as 
the pseudo-SC/ (3) do not have this property. In their re- 
cent investi gati on of the onset of rotational motion, Zuker 
et. al. l26l] introduced a preliminary formulation of 
an approximate quasi-S(7(3) symmetry. It would be in- 
teresting to study whether the quasi-S'[/(3) contains a 
similar property. 



Fig.l presents a highly non-trivial result because the 
PSM, as described above, is not built on any explicit 
SU(3) symmetry, and no free parameters have been ad- 
justed to obtain such a symmetry. Nevertheless, the spec- 
tra, the electromagnetic transition rates, and the wave 
functions of the PSM agree nearly perfectly with the 
FDSM SU (3) results, from the grou nd-band to states of 
high spins and high excitations This strongly sug- 
gests that the projected deformed-BCS vacuum in the 
PSM could be at the microscopic level close to the S- 
D core in the FDSM, provided that the S-D pairs must 
be redefined in a multi-major-shell space, since the one- 
major-shell is not large enough to accommodate so many 
active nucleons (here n° = 32 and n^. = 16). 



In this regard, there is a conceptual distinction be- 
tween the SU(3) symmetry in the one-major-shell FDSM 
and that emerged from the PSM. In the former, the 
SU(3) symmetry arises entirely from the normal-parity 
nucleons, and the abnormal parity orbit enters only im- 
plicitly through the Pauli effect and the renormalization 
of the parameters. In the latter, the SU(3) symmetry 
arises from the explicit dynamical participation of both 
normal and abnormal-parity nucleons of many shells. Ig- 
noring the direct contribution from the abnormal-parity 
nucleons in the one-major-shell FDSM is a sacrifice for 
having an exact symmetry (which forms a closed Lie al- 
gebra). In practice, this turns out to be not a bad ap- 
proximation for collective motions because in a single ma- 
jor shell, there is only one single- j level with abnormal- 
parity, which does not have much quadrupole collectiv- 
ity compared to the normal-parity contributions |27| . 
When going to a multi-major-shell space, the situation 
will change. A bunch of abnormal-parity levels, which 
are located just below the normal-parity levels, will open 
up. This means that, when one redefines the FDSM-type 
coherent S-D pairs in a multi-major-shell space, the con- 
ceptual distinction of the SU (3) symmetry between the 
FDSM and the PSM will be eliminated. 
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B. Collective D-Pair Excitations 

To see further the relationship between the two mod- 
els, let us ignore the terms with r ^ 2 in the most general 
FDSM Hamiltonian H12|l . This is a reasonable approach 
because from the multipole expansion point of view, there 
is no multipole interactions with r — odd without consid- 
ering the parity admixture. One may notice that there 
is the monopole-monopole interaction (r = 0) in the 
FDSM, which is not included in the PSM. However, it is 
well-known that this interaction only affects the nuclear 
total binding energy but does not have much influence 
on the excitations. On the other hand, it can be added 
to the PSM if necessary. Thus, we have for both models 
the interactions of the monopole- and quadrupole-pairing 
plus the quadrupole-quadrupole type. 

The PSM Hamiltonian contains operators written in 
the ordinary shell model basis, whereas the FDSM ones 
in the k-i basis. In order to compare them, we transform 
the FDSM operators in Eq. l|51 ITT|l back to the ordinary 
shell model basis by rewriting them in terms of c£, and 
c Q . We find that 



P 



2/i 



V {2) c f c , 



^ — q c l c a 



(14) 



J (2) JJ, 

2 £ ^ /j. a a a 

ol.ol' 



with 



(2) 



= \/2n kl c 211 , ,(-) 



i k j 
i k j' 
2 2 



Comparing Eq. (14) with Eq. @, we see that the defini- 
tions of St, Dt and P 2m in the FDSM are similar to that 

of P>, P^ and Q M in the PSM, except that the model 
spaces for the two models are different. However, within 
one-major-shell, the operator S' (S) is exactly the same 
as Pt (p) in the PSM. For the £>-pair and P 2 opera- 

(2) 

tors, although the coefficients Q^a' different from 
Qnaa' appearing in the corresponding operators and 

in the PSM, their physical meanings are the same. 
The FDSM Hamiltonian may thus be considered as an 
one-major-shell version of the PSM Hamiltonian, with 
the approximations of P^ and being replaced by DJ^ 
and P 2/J and the s.p. energy splitting being ignored. In 
other words, if, in the multi-major-shell case, the sym- 
metry constraint is released from the FDSM and the s. p. 
energy splitting is considered, Sr , -DI and P 2 ^ in the 
FDSM should return back to the version of P\ P,t and 

Qfj, operators in the PSM. So is for the Hamiltonian. 

From the above analysis we may conclude that the 
PSM and the FDSM are just two approaches to solve 
an effective Hamiltonian of a common form. In order 



for them to be applicable to heavy nuclei, approxima- 
tions have to be made in each model. The PSM is based 
on a multi-major-shell space so that it can describe the 
nuclear rotational motion microscopically through a dy- 
namic participation of many particles. However, it can 
afford to do so only for a truncated configuration space 
that includes only a BCS vacuum plus a few qp excita- 
tions. As a sacrifice, this truncation do not include the 
collective modes such as (3- and 7- vibrations. In con- 
trast, the FDSM aims at nuclear low-lying collective ex- 
citations. It can afford to do so only when the model 
space is reduced to one-major shell like in the conven- 
tional shell model. But that is still not enough for heavy 
nuclei. Additional approximations to further reduce the 
configuration space down to the symmetry detected S-D 
subspace are necessary. These approximations exclude 
the s.p. excitations. 



IV. CONSTRUCTION OF HEAVY SHELL 
MODEL 

Having realized that the projected deformed-BCS 
states exhibit the SU(3) symmetry, and that the 
collective excitations may be approached by the im- 
pair excitations, with the D-pair operator defined as 
the quadrupole-pair operator P^ in a multi-major-shell 
model space, we propose a multi-shell shell model: the 
Heavy Shell Model. The essence of this proposal is to 
adopt the truncation scheme for the collective modes, 
which was discovered by the FDSM, into the PSM to en- 
rich the shell model basis. This essentially combines the 
advantages of both models. 

In fact, to incorporate both s. p. and collective exci- 
tations, there are in principle two alternatives: One can 
either extend the FDSM by adding the PSM qp trun- 
cation scheme on top of the FDSM collective states, or 
extend the PSM by including the D-pair collective ex- 
citations into the PSM vacuum. However, the FDSM is 
a severely truncated one-major-shell shell model dictated 
by symmetry, and in this sense, it is not as microscopic as 
the PSM. Although it is quite successful in the descrip- 
tion of low-lying collective motions, the FDSM is just an 
effective theory. For practical applications, large renor- 
malization effects must be embedded in the parameters 
of the FDSM Hamiltonian, which have to be determined 
phenomenologically. Therefore, we choose to construct 
the Heavy Shell Model based on the extension of the 
PSM. 

The main ingredients of the Heavy Shell Model are as 
follows: 

1. Keep the multi-major-shell basis of the PSM as 
the model space, using the PSM pt and P 2 ' op- 
erators to describe the coherent S- and Z?-pairs in 
the multi-major-shell configurations, and construct 
the intrinsic collective excitation states by P 2 act- 
ing on the deformed-BCS vacuum; 
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2. Carry out the shell model truncation by selecting 
a few single-qp states near the Fermi surfaces plus 
a few D-pair excitations, and perform angular mo- 
mentum and particle number projection to obtain 
a shell-model basis in the laboratory frame; 

3. Keep the PSM Hamiltonian to be the effec- 
tive Hamiltonian, but allow adding more multi- 
pole interactions and/or readjusting the interaction 
strengths if necessary; 

4. Utilize the algorithms developed in the PSM to 
carry out calculations for all the necessary matrix 
elements, and diagonalize the Hamiltonian in the 
truncated shell model space. 

Let us now discuss each of the items in more detail. 



A. The Basis States 

We have demonstrated that the projected deformed- 
BCS vacuum |$) in the PSM is nearly identical to 
the FDSM-5?7(3) intrinsic ground state (n eff , 0) irrep. 
Furthermore, we have indicated that the FDSM S- 
and D-pair and the P^^ operator are, respectively, 
the symmetry-constraint one-major-shell version of the 
ordinary monopole pair P, quadrupole pair P^, and 
quadrupole operator in the PSM. It is therefore nat- 
ural to believe that in a multi-major-shcll space, the 
FDSM S- and D-pairs are nothing but the P- and im- 
pairs if we abandon the symmetry requirement. So is 
for the quadrupole operator. In the FDSM, all known 
types of low-lying collective excitation can be obtained 
by acting Z)t on the FDSM-S'[/(3) intrinsic ground state, 

(St + PDl) N \0) (see Table I). Combining these facts, 
the collective excitations of the HSM in a multi-major- 
shell space may be constructed by replacing, respectively, 
(Si+0Dl) N \O) and £>t m Ta bl e I with |$) and pt denned 
in the PSM. 

Hereafter, we will continue to use the FDSM notations 
S and D for the pair operators. One should bear in mind 
that they have been redefined as 

^ - ^£44 

a 

D l = ^ = £ 44,. (15) 

a, a' 

The intrinsic collective states can then be expressed as 
|* c ) = \N D (np,n^K)) =IIl>jJ*) 

i=0 

= (Dl) n ?(DtDl 2 )^(Dl)%\<Z>). (16) 

Eq. (|16J) provides the microscopic meaning of the quan- 
tum numbers np, n 7 and k in a very clear manner: the 



phonon number appearing in phenomenological models 
is nothing but the total number of D-pairs 

N D = n p + n 7 + k/2. (17) 

The basis of the HSM can be constructed by adding 
qp-excitations on top of the collective intrinsic states; 
the formalism is the same as that used to build the 
PSM bases in Eq. (1), but the simple BCS vacuum |$) in 
Eq. (1) is now replaced by a more correlated one, |3> c ). 
The general expression of the HSM basis in the labora- 
tory system can be written as 

\qcIM) = P^l^c) 

\%c) ee nn ( is ) 

where n v q (n*) is the qp number of neutrons (protons), 
and the indices q and c stand for the qp and the collective 
vibrational configurations, respectively. 

The HSM basis (|18fl contains both s. p. (qp excita- 
tions) and collective (D-pair excitations) degrees of free- 
dom; the deficiency of lack of collective degrees of free- 
dom in the original PSM is redeemed. Moreover, the 
basis i)18fl is expected to work also for the transitional 
(or weakly deformed) nuclei, which is beyond the origi- 
nal PSM territory. This is expected because it is known 
from the FDSM that the collective states of transitional 
nuclei can be described as a mixture of different SU (3) ir- 
reps The original PSM uses only the deformed-BCS 
vacuum (the ground state of the SU(3) irreps), and thus 
does not contain such a mixing mechanism. This is why 
the PSM becomes less and less valid when going away 
from the well-deformed region. The HSM basis (|18fl now 
contains all possible 517(3) irreps, since its labels {rip, 
n 7 , k) are in one-to-one correspondence to that of the 
SU{3) irreps (A,/x) 

For spherical nuclei, the basis should be constructed 
separately since in the spherical case, the rotational sym- 
metry is restored so that no distinction can be made be- 
tween the intrinsic and the laboratory system. In the 
spherical limit, the BCS vacuum is just the 7 = ground 
state, which, in the FDSM, corresponds to the intrinsic 
state with (3 — ► (see Table II). Let us denote this BCS 
vacuum as |$o) to distinguish it from the deformed one. 
The deformed mean- field (e. g. the Nilsson scheme) in 
the spherical case is also reduced to the spherical shell 
model level scheme so that a qp state can be labeled 
by {N£jm}. Using {jm} for short, a qp state can be 
expressed as ajtJ$o) (a — v,ir), which is no longer a su- 
perposition of all possible angular momentum states but 
has a definite angular momentum. Likewise, a T>-pair ex- 
cited state 7)t |<J> ) is now a 2 + state. Therefore, there is 
no need for performing angular-momentum projection. 
The basis with a given total angular-momentum 7 can 
be directly obtained by the angular-momentum coupling. 
The general expression of the HSM basis in the spherical 
case can be written as follows: 
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\qcIM) = P N \^) 

= [a^)®^^)]^) (19) 



with 



4 c uK) =\AUn»)®AUnl) 



M„ 



(20) 



4 f k)= i n a i\ m«=y. c m n 4m, ( 2i ) 

i=l [m] i=l 

VkMS N °)= i II =E C M (RcM c ) J] £>t fc (22) 



fc=i M 



fc=l 



where ^ JaMa (n,) and Djj M (A^) are the creation 

operators of qp's and Ap D-pairs coupled to 
angular-momenta JgM° and R C M C , respectively. The 
short-hand notation [m] represents the configuration 
{mi,TO2,- • • ,m n °} of the qp's. Similarly, [/i] is for the 
_D-pairs. 

The amplitudes C{ m }{J°M°) in Eq. (|2*T)l can be easily 
calculated by using the standard shell model technique. 
As a matter of fact, all spherical nuclei lie very close 
to the doubly closed shell. Therefore, only a few qp's 
from the j-orbits around the Fermi surfaces need to be 
considered. The amplitudes Cnj (R c M c ) in Eq. 1221) are 
difficult to obtain in this way because multi-major shells 
are involved due to the Z?-pair collectivity. However, they 
can be easily evaluated through the d-boson CFP's (Co- 
efficients of Fractional Parentage) if the Z?-pairs are re- 
garded as bosons. It should be noted that this is not a 
boson approximation since the Dt's in Eq. JUJ remain 
to be fermionic operators. It is just a symmetry detected 
truncation to select only those states that are symmetric 
with respect to interchanging any Z?-pairs. 

A comparison between the deformed basis and the 
spherical basis is shown in Table II. It is interesting 
to see that the D-pair excited states, which, in the 
deformed case, produce the (3-j bands after angular- 
momentum projection, degenerate in the spherical case 
into a spherical vibration-spectrum. The spherical 
vibration-spectrum has much less independent states 
than the rotational spectrum. It is so simply because 
when P — > 0, most of the states obtained from the 
angular-momentum projection are not linearly indepen- 
dent due to the spherical symmetry nature. They are 
highly over-complete. For instance, the f3- and 7-bands 
reduce to a 2 + state when f3 — ► 0. This is why we should 
construct the spherical shell-model basis differently. It 
can be seen that even without configuration mixing, this 
collective basis (as shown in Table II) can already give 
the essential features of the low-lying collective modes 
found in the spherical nuclei. We therefore feel confident 
that the HSM basis is going to work well. 



Of course, having built a shell-model basis is not 
enough. The central theme that follows is whether it 
is possible to find an efficient truncation scheme within 
the constructed basis so that the calculations for heavy 
nuclei become feasible. 



B. Basis Truncation 

To estimate the size of the deformed intrinsic basis 
(|18[) . let us denote the maximum qp number, the num- 
ber of selected s. p. states, and the maximum number 
of the excited D-pairs as n%, n° (a = v,tt), and N m , 
respectively. The total dimension is the product of the 
dimension of neutron qp states, that of proton qp states, 
and that of the D pairs 



Da 



A^KX) x 2?$«X) x D$(N m ). (23) 



im V q 5 si im \ q ' s , 

For each of these three terms, it can be shown that 



(E„ is in step of 2) 



D 



N m N m -N D 
{Nm) = E E (■ 



Nf? +4A^+3+<5 



N"=0 iVK=0 



, Viy5 a +4iV5+3+^ 



where D)^} (nZ, nZ), with a = v or 7r, are the dimensions 



(24) 



im v' "q ■ 

of the qp states for neutrons or protons, D y ^ l {N m ) is the 
dimension of the intrinsic collective states, with S a — 
1 when Nf) = even and 8" = when Ag = odd. In 

Eq. I|24|l , D^ (N m ) is for separated neutron and proton 
vacuum [15J . If we treat the neutron and proton vacuum 
as a single coupled BCS vacuum, the double-summation 
in Eq. is reduced to a single one and this can be 
easily summed up: 

D^(N m ) = i[ (A m + l)(A m +3)(2A m +7)+3<5 Nm ] , 

(25) 

with <5jv m =1 or depending on N m = even or odd. 

According to the experience of the PSM calculations, 
it is sufficient in most cases to take nZ < 2 and nZ < 4 for 
both neutrons and protons. One can then obtain that the 



dimension of the qp states D^ n {n v q , n v s ) x £)^(n£, nj) < 
49. If Eq. (|25|l is used for estimating dimension of the 
collective states, D^ m (N m ) changes from 1, 3, 7, 13 to 161 
when the maximum excited _D-pair number N m varies 
from 0, 1, 2, 3, to 10. Thus, the total dimension 
Am ranges from 49, 147, 343, 637, • • • , to 7889. 

For well-deformed nuclei, the Z?-pair excitation energy 
is about 1 MeV (i. e. the band head of the first f3-j 
band). States with more D-pairs are excited higher in 
energy. Furthermore, well-deformed nuclei are very close 
to an SU(3) rotor. States with different Ap, which cor- 
responds to different 5C/(3) irreps, hardly mix with each 
other although the mixing may not be zero since the sym- 
metry is not perfect. Thus, for the dominant low-lying 
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TABLE II: Comparison of collective basis for deformed and spherical Nuclei. * 



Deformed Nuclei 




Spherical Nuclei( (3^0) 




FDSM SU(3) 


HSM 


Rotational 


FDSM 


HSM 


Vbrational 


Intrinsic States 


Intrinsic Bass 


Bands 


SU(2)States 


Basis 


Spectnm 


(s''+ poj)" |o) 


I*) 


Ground Band 


(sY |o) 


Icb \ 
1*0/ 


+ 0?ko 


d:( S >+VdI) n -'\o) 


D i*> 


P- Band 
y- Band 


d;(sT|o) 


oil*,,) 


2 + 1 ftco 


u o u o ^ + P u o > l u / 




(020) Band 


V\ u ) JoV^ ) | u / 


> Jo 1 0/ 


+ 


d^dUs\ bz>V" 2 Io> 
/>.'/> (s • Pd;)"- 2 |o) 


d'd'\$>) 


(Oil ) Band 
t^uuz ) Ddno 






2 + 2/ko 


l> l> . •: S • PD* j"" 2 |0) 


D'p'\<t>) 


(200) Band 


[(D%(.S f f |0) 


tpVitl*,) 


■1 + 

4 


D tD o tfl o t (* t + K)" 3 |°) 


d'd'p')®} 


(030) Band 


urffivT |0) 


[p f )']° 1*) 


o + 


D 2 t DjDj(s , + P£> T ) W 3 |0> 
DlD\Dl (S f + PD *)' v ^|0) 




(021) Band 
(012) Band 


i(rf)%(sT |o) 


[P f )tl*») 


2 + 




d'd*d'I<[>) 


(003) Band 


[fB t )t(S , f |0> 




3 + 3/zco 


D 2 f D: 2 Do(S f + PDoV^lo) 


dTj'/j'J*) 


(210) Band 




[p + ) 3 ];i*.) 


4 + 


£>*£>* D*,(sV PZ)*)*" 3 |0) 




(201) Band 


Krtff^f |0) 


[p , y , ];i*> 


6 + 



* The first column lists the FDSM SU (3) intrinsic states. The second column is the HSM collective basis for the deformed 
nuclei, which is obtained by replacing (S^ + {3Dq) n \0) with the BCS vacuum 1$}. After projection, each intrinsic state 
produces a rotational band labeled as /3, 7, and (n 7 , up, k/2) listed in the third column. The FDSM 5(7(2) states and the 
HSM collective basis for spherical nuclei are listed in the 4th and 5th column, respectively. The last column lists their 
corresponding spin-parity and excitation energies relative to the BCS vacuum. Note that the D-pairs in the SU(2) states 
have been modified to commute with 



collective states the number of the excited ZJ-pairs, Nd, 
should not be large. We expect that taking N m — 2 to 4 
is already good enough to produce all the low-lying col- 
lective vibrational states. Thus, the total dimension D im 
is of the order of 10 2 to 10 3 for well-deformed nuclei. 

For weakly-deformed nuclei, mixing between different 
5(7(3) irreps must be stronger, and therefore, it requires 
a larger N rn to account for such a mixing. This can in- 
crease the basis dimension drastically. However, transi- 
tion from the deformed to the spherical region behaves 
like a phase transition [2^, |3(j , which happens suddenly 
within a small interval of nucleon number. We therefore 
expect that as we pass through the weakly-deformed re- 
gion, nuclei may quickly enter into the spherical region 
before the N m number becomes too big (e.g. N m > 10). 
Thus, calculations for weakly-deformed nuclei, though 
harder, are tractable. 

When including the scissors mode in the calculation, 
we can use Eq. (|24[) for estimating the basis size. The 
maximum dimension in this case would rise by one to two 
orders of magnitude. Namely, the dimension would be- 



come 49, 294, 1127, 2646, • • • , 2xl0 5 for N m = 0, 1, 2, 3, 
• • • , 10. Moreover, the dimension should be multiplied by 
a factor of /„ p , since, after the neutron and proton states 
have been built with good angular momentum, there are 
many ways of coupling these neutron and proton states 
to a given total angular momentum /. It can be shown 
that for a given total angular momentum / < 2I cut , 

fn P (I) = -(4/cut + 1 - 37) + {I cut + 1), (26) 

where I cut is the cut-off angular momentum (i. e. the 
highest angular momentum that are considered for neu- 
trons and protons). This factor is of an order of 10 2 for 
Icut ~ 10. Thus, the total dimension Di m would be of the 
order of 10 5 to 10 6 for N m — 3 to 4. However, given the 
fact that different collective modes usually have different 
energy scale and different symmetries, there is still room 
for a further reduction in dimension. It is known that 
the excitation energy of the low-lying scissors mode in the 
well deformed region is about 2 MeV higher than the low- 
lying /3~7 bands. Furthermore, as shown in Refs. [l4lll5| . 
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the scissors mode is physically caused by the relative mo- 
tion of separated n- and p- vacuum (with N rn — 0), not by 
the D-pair excitation. Therefore, the interplay between 
these different collective modes should be small and may 
thus be studied separately. One can thus set N m = 
when studying the scissors vibrations. If we include qp 
excitation states up to 4-qp and no D-pair excitations in 
the basis, the dimension can be reduced to about 5xl0 3 . 

For spherical nuclei, the dimension should be estimated 
differently. What Eqs. (%3|l and (|2*5)l show us is the total 
m-scheme dimension, Di m , in the intrinsic frame. What 
we are actually interested in is the dimension Di m (I) 
at a given angular momentum /. In the deformed case, 
Di m (I) is equal to D im because after projection, every 
intrinsic state can in principle contribute one state to an 
I. In the spherical case, D im (I) is smaller than D im since 
the basis is now constructed by angular-momentum cou- 
pling in the laboratory frame. Eq. (|23|l may be utilized 
to estimate the m-scheme dimension, but the formula for 
counting the collective basis states in the spherical case 
should be replaced by 

D&(N m ) = £ ( Nd + 4 ). (27) 

N D =0 \ ' 

This formula is obtained because the collective basis 
states are constructed differently for the spherical case, 
and we have imposed a constraint that only those ba- 
sis states that are symmetric with respect to the inter- 
change of D-pairs are selected. According to Eq. (12711 . 
for N m = 0,1,2,3,- •• , 10, the dimension D^(N m ) is 
respectively 1, 6, 21, 56, • • • , 3003, which are consider- 
ably larger than the dimensions in the deformed case. 

An estimation for the maximum dimension of Di m (I) , 
denoted as D™£ x , can be made by using the relation 
A m = Eo max D im {I)(2I + 1): 

T\max p r) ^ 4£>j m 

^irn ~ niU im - (7max+2) 2 

T = T u + V 4- T c 

^max ^max 1 max 1 maxi 

where I max (u = v, 7r, c) are the maximum angular mo- 
mentum that the neutron qp's, the proton qp's, and the 
excited D-pairs can reach. Thus, knowing the m-scheme 
dimension Di m one can estimate the maximum dimen- 
sion D™ x . The reduction factor i? 7 = [2/{I max + 2)] 2 
makes D r {!£ x two to three orders of magnitude smaller 
than the m-scheme Di m . 

On the surface it seems that the basis for the spheri- 
cal case is smaller because of the reduction factor Ri in 
Eq. H28fl . it is in fact not true. The problem lies in the 
fact that in the spherical case, the density of s. p. states 
around the Fermi surfaces is much larger than that in 
the deformed case since each j has a 2 j + 1 degeneracy. 
The number of qp states rig must be around 30 instead 
of 4 adopted in the deformed case. This increases the m- 
scheme dimension on the qp sector enormously. Keeping 
n q < 2, the dimension of qp basis will increase from 49 



in the deformed case (n a s = 4) to an order of 10 5 , leav- 
ing almost no room for collective basis. This is another 
reason why we can not continue to use m-scheme basis 
for the spherical case. Using /-scheme only D™£ x should 
be concerned, which greatly reduces the dimension by a 
factor of Rj. The total basis dimension varies from 1125, 
5879, 17742, 4xl0 4 , • • • , to lxlO 6 for N m = 0, 1, 2, 3, • • • , 
to 10, when I^ ax = l v max = 12 is assumed. 

It is known that the low-lying spherical vibrational en- 
ergy hw is about 0.5 MeV, which is roughly the Z?-pair 
excitation energy in spherical nuclei. The s. p. excitations 
start to have influence at about 2ftw and become signifi- 
cant around Sfku, where the characteristics of vibrational 
spectrum are strongly disturbed. Thus on the collective 
sector, taking N m — 3 (3fku excitation energy) may be 
a reasonable choice. If the number of qp's is kept to be 
n a q — 2 (a — v, 7r) to give the maximum 4-qp states, the 
dimension one has to deal with is of the order of 10 4 . 
In fact, spherical nuclei lie often near the doubly closed 
shells and have large binding energies. The 4-qp states 
may already be rather high and have only little influence 
to the low-lying states. If we ignore the 4-qp states, then 
the maximum shell-model dimension can be reduced to 
488 for N m — 3. Even for N m = 4 the maximum dimen- 
sion is only 908. Such a basis is clearly manageable. 

C. The Effective Interactions 

The HSM Hamiltonian can be generally written in the 
following form 

H= £ H a + H„ w , H vw = -J2 xl X J Q'i ■ Ql 

a — v, 7r A 7^ 1 

A jt 1 

(29) 

This rotational invariant Hamiltonian is the same as that 
in the PSM || (see also Eq. iJBJ), except that the multi- 
pole interactions are now extended to include not only 
the quadrupole but also monopole, octupole, and hex- 
adecupole terms. All these operators are defined in the 
multi-major-shell space: 

= Q^iaa' c a c a'i 

a, a' 

Q%> = {a\r x Y Xll {e^)\a') , (30) 
A = 0,2,3,4. 

The octupole and hexadecupole interactions have been 
employed in their Hamiltonian by Chen and Gao [3ll | in 
dealing with the actinide nuclei with projection. This 
type of schematic interactions l|29l) works for the struc- 
ture calculations surprisingly well despite its simplicity. 
In fact, it has been shown by Dufour and Zuker [32] that 
these interactions simulate the essence of the most im- 
portant correlations in nuclei, so that even the realistic 
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force has to contain at least these basic components im- 
plicitly in order to work successfully in the structure cal- 
culations. Therefore, we find no compelling reasons for 
not using these simple interactions. Of course, the model 
is open to adoption of any realistic forces. 

The s. p. energy term Hq can be simply taken from 
the Nilsson scheme at zero deformation although other 
schemes such as the Woods-Saxon may also be adopted 
if there is an advantage. The parameterization in the 
Nilsson scheme has been well established for the /3-stable 
mass regions (see, for example, Ref. jsH)- F° r the ex- 
otic mass regions, the standard Nilsson parameters may 
not be valid and improvement may be necessary |34l l35| . 
In this regard, new experimental data that can provide 
information on the single-particle energies in the exotic 
mass regions are very much desired. The Relativistic 
Mean Field Theory (RMF) [H H3 could also be help- 
ful in providing a s. p. energy scheme for the exotic mass 
regions where no data are available for determining a phc- 
nomenological mean field. The RMF may have a better 
extrapolation power than other phenomenological mean 
fields because a single set of parameters of the RMF is 
able to fit nuclear ground state properties from light to 
heavy nuclei reasonably well. In practice, it may be a 
convenient approach to use the RMF s. p. energies as a 
reference to adjust the Nilsson parameters in the HSM. 

The parameters (A = 2,3,4) in Eq. J2U) are 

determined by the self-consistent relation with defor- 
mation parameters, in the same way as in the PSM 
0,1^. The monopole-monopole interactions (A=0) in- 
cludes three terms: xl n " ( CT = v > 7r ) an< ^ X^} n u n -K- 
It has been shown that the monopole-monopole inter- 
actions are important for the nuclear mass calculations, 
in particular for those nuclei lying far from the stabil- 
ity line |38l l3Si | . They can be viewed as an average way 
of accounting for the N-Z dependence in s.p. energies. 
However, for a given nucleus, the monopole-monopole 
interactions contribute a constant energy only. Thus, if 
we are not interested in calculating the absolute energies, 
these terms can be ignored. The strengths of monopolc- 
pairing, Gjjj, and quadrupole-pairing, Gq, depend on the 
size of single-particle space. They are inversely propor- 
tional to the mass number A. 

For the well-deformed mass regions where calculations 
have been extensively performed by the PSM, the inter- 
action strengths proposed for the HSM should be simi- 
lar if the same size of single-particle space is employed. 
However, when the HSM is applied to other mass re- 
gions such as the transitional or spherical region, these 
strengths may need to be readjusted. In particular, the 
self-consistent relation used for determining the multi- 
pole interaction strengths will break down when the ba- 
sis deformation becomes zero. The strengths for these 
cases have to be studied separately. After all, Eq. lt2D|) 
is an effective Hamiltonian in nature, and is subject to 
vary when the model is applied to different mass regions. 
Nevertheless, a smooth variation in parameterization is 
expected since the model space is sufficiently large. 



D. Evaluation of Matrix Elements 

Having had a tractable basis and a reliable effective 
Hamiltonian, the remaining task is to diagonalize the 
Hamiltonian in the basis to get the eigen-energies and 
eigen-functions, and then to use the obtained wave func- 
tions to calculate the observables. To do so, one must 
know how to evaluate the matrix elements in the pro- 
jected basis. The projection techniques have already 
been well developed by the PSM based on the pioneer- 
ing work of Hara and Iwasaki |40j . The extensive discus- 
sion about the details of the projection techniques can be 
found in the PSM review article and the PSM com- 
puter code |41| . Here we emphasize only on how to deal 
with the new ingredient: the D-pairs and the spherical 
basis. 

Let us first discuss the deformed case. Suppose that a 
one-body tenser operator of A-rank, T^, and the eigen- 
functions |^^f) are given in the laboratory frame: 

= E T ^a'cic a , , ) = E F^Pj&l^c) (31) 

v.v' qc 

where T* aa , = (a\T\fj_\a') is known. \Q qc ) is the intrin- 
sic basis defined in Eq. lfT5|l . and the amplitudes Fl c are 
obtained by solving the eigen- value equation, Eq. 10} . Al- 
though we take the tenser operator in Eq. (|31|l as exam- 
ple, the following discussion applies equally well to the 
pairing operators if we change cj,,c a < in Eq. I|31l) to c i a c i a i 
or c a c a > and rewrite T^ aa , = (Q\f^\aa') or (aa'|T AAl |0). 
Sandwiched by the wave functions, the matrix element 
of T\ M in the laboratory frame can be expressed as 

(<,|T V K,) = 

E ^>,c(^'o'|pSJa/SI^c). (32) 

q'c'qc 

By applying [NT^] = ANT\^ where AN = (±2) for 
the electromagnetic multipole (pairing) operator, it can 
be shown that 

r>l'N' rf, f,IN _ cN' fil'M' 

x E ^iK'-u' ^ii'TXji'PjV-n'K- (33) 

Inserting Eq. O into Eq. (|52"jl . we obtain 

1^1*5,) = C^'^'^xW), (34) 

with 

(* / UTaII* 7 ) = 5jv + AiyE F l' c! F gc C IK'-u',\u' 
fi' q'c'qc 

x (^M^Pk^K^) (35) 

and 

{fq'AT^Pk^K^qc) = 2 -^- J d<t>dn e iN <f> 

x^ir(fi) (^vlTVM^I^c). (36) 
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Here, R{fl, 4>) = R{VL)e~ l ^ N is a rotation operator in the 
4-dimensional (coordinate plus particle-number) space. 
The problem is then reduced to evaluate the matrix el- 
ements having the general form (& q ' c > \T\ V R(Q,, <f>)\^ qc ). 
These are the matrix elements for a one-body tenser op- 
erator between the intrinsic state ($ g ' c '|, constructed in 
(fTHft . and the rotated one R(S},<p)\$ qc ). 

In order to evaluate the matrix elements, it is conve- 
nient to transform the operators cj, (c a ), contained in 
Txu, and in the D-pairs (which are embedded in the basis 
|$ gc )), into the qp representation {a\a}. The transfor- 
mation matrices (of the Hartree-Fock-Bogoliubov type) 
are known: 

with XJ av = W av u v , V av = W av v u , (3' 

where u and v are the occupation amplitudes in the BCS, 
and PFqi/'s the Nilsson wave functions of the deformed 
s. p. level v. 

The evaluation of ($ q , c ,\f Xtl 'R(n, 4>)\<& qc ) in Eq. ^ is 
equivalent to evaluate the following general expectation 
quantity with respect to the BCS vacuum state |$): 

($ 9 v \T\fi'R(Cl, </>)\§ qe ) = ($|--- D^'D^y ■ ■ ■ a v2 ,a vV 
x f x ^R(n, 4>) a^aU ■ ■ ■ D^D+a ■••]$). (38) 

After applying the transformation from the spherical to 
the deformed-BCS basis (Eq. (|37|l ). The operators ap- 
pearing in the above expression, such as T\ M and D^, are 
generally linear combinations of products of the qp oper- 
ators and a. The problem eventually becomes the one 
of performing contractions for a series of qp creation and 
annihilation operators 

($| • • • a D2 /a D1 / ■ • • a v2 'a v v • • • a\, 2 a T1 R(VL, (j>) 

x 4i4 2 • • • 4iaL • ••!$)• (39) 

In performing the contraction calculations, a generalized 
Wick-theorem is used. The techniques of carrying out the 
contractions with the rotation operator R are available. 
It has been shown 0,^3 that the problem can be reduced 
th evaluation the following three basic elements 

A vv < = (a v [Q]al,), B vvl = (a„a„/[fi]>, C vv > = ([n]a^a^,), 

(40) 

with [f2] = R/(R), and (• • •) is the short-hand notation 
of vaccum expectation, for example, (O) = ($101$). The 
contraction can then be calculated by the following the- 
orem 

(a n > ■ ■ -ai/[0]a} • • ■ a^) = 

[f] 

^(±)(B) fe ~ (n ~™' )/2 (CT~ 2fc (A) fc (41) 

k—ra p 

(n + n = even, m — max{0, — - — }) 



where J2 P 1S a "permuted sum" with all possible com- 
binations of pairs of indices, and ± is the permutation 
parity. Details for these calculations can be found in the 
PSM review article @. 

For the spherical case, the matrix elements can be eval- 
uated in the same way. The only difference is that there 
is no need for doing angular-momentum projection since 
the basis \& q c ) defined in Eq. I|19|) is constructed in the 
laboratory frame. Therefore, in the spherical case, the 
projection operator P\f K should be replaced by P N . As 
a consequence, Eqs. (13511 and (|36|) become 

pi,' q' c'qc 

x {^'\f Xtl ,P N \^ q f) (42) 

and 

x<^: c r|f v e-^|^f). (43) 

After transforming and the D-pairs embedded in 
|$g/^f ) and |$q^) into the qp basis, the evaluation of 

matrix element \fx tl >e~ lN ' l '\§ I q c I ) is again reduced 

to the contraction calculations like the type of Eq. ll3"?)l) , 
except that in this case, R is replaced by e~ lN ^ . 

V. SUMMARY 

In this paper, a multi-shell shell model for heavy nuclei 
is proposed. For performing a shell model diagonaliza- 
tion involving several major shells in the model space, 
we seek an efficient truncation scheme. The new Heavy 
Shell Model can be viewed as an integration of two ex- 
isting models: the Project Shell Model and the Fermion 
Dynamical Symmetry Model. The PSM is an efficient 
method for the high-spin description of rotational states 
built upon qp-excitations, but it is not a practical method 
for the low-spin collective vibrations. In contrast, the 
FDSM provides a well-defined truncation scheme for all 
known types of low- lying collective vibrations, workable 
from the spherical to the well-deformed region, but it is 
lack of the necessary degrees of freedom of single particle 
excitations. The idea proposed in the present paper is 
to combine the advantages of both models. To construct 
the shell model basis, we follow the FDSM discovery that 
the intrinsic collective states can be built by applying im- 
pairs onto the BCS vacuum state, and employ the PSM 
qp truncation scheme combined with the projection tech- 
niques. In this sense the model goes beyond the tradi- 
tional one-major-shell shell model, yet the calculation is 
tractable for heavy, and even for superheavy nuclei. 

Given the past success of the PSM and the FDSM 
in their own applicable regimes, which has been docu- 
mented in the literature, we expect that the new model 
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can work reasonably well for heavy nuclear systems where 
traditional shell-model calculations are not feasible. This 
model should be capable of describing the low-excitation 
collective vibrations of all known types, including frag- 
mentations due to the quasi-particle mixing. It should be 
capable of applying to the high-spin regions where quasi- 
particle alignments play an important role. It should also 
be capable of treating weakly-deformed nuclei across the 
transitional to the spherical region. We conclude that the 
development of the Heavy Shell Model may open possi- 
bilities of shell model calculations for heavy nuclei to a 



much wider range of nuclear structure problems. 



Acknowledgments 

Research at the University of Notre Dame is supported 
by the NSF under contract number PHY-0140324. Re- 
search at Physics Division of National Center for Theo- 
retical Sciences is supported by National Science Council 
in Taiwan. 



[1] P. Ring and P. Schuck, The Nuclear Many-Body Problem 
(Springer, New York, 1980). 

[2] S. Frauendorf, Rev. Mod. Phys. 73, 463 (2001). 

[3] K. Langanke and G. Martinez-Pinedo, Nucl. Phys. A 
673, 481 (2000). 

[4] E. Caurier, A. P. Zuker, A. Poves and G. Martinez- 
Pinedo, Phys. Rev. C 50, 225 (1994). 

[5] CA. Ur, D. Bucurescu, S.M. Lenzi, G. Martinez-Pinedo, 
DR. Napoli, D. Bazzacco, F. Brandolini, D.M. Brink, 
J. A. Cameron, E. Caurier, G. de Angelis, M. De Poli, A. 
Gadea, S. Lunardi, N. Marginean, MA. Nagarajan, P. 
Pavan, C. Rossi Alvarez, and C.E. Svensson Phys. Rev. 
C 58, 3163 (1998). 

[6] K.W. Schmid and F. Griimmer, Rep. Prog. Phys. 50, 731 
(1987). 

[7] M. Honma, T. Mizusaki, T. Otsuka, Phys. Rev. Lett. 75, 

1284 (1995). 
[8] T. Otsuka, Nucl. Phys. A 704, 21c (2002). 
[9] K. Hara and Y. Sun, Int. J. Mod. Phys. E 4, 637 (1995). 
[10] C.-L. Wu, D. H. Feng and M. Guidry, Adv. Nucl. Phys. 
21, 227 (1994), ed. J.W.Negele and E. Vogt (Plenum 
Press, New York-London). 
[11] E. Caurier, J.L. Egido, G. Martinez-Pinedo, A. Poves, 
J. Retamosa, L.M. Robledo and A. P. Zuker, Phys. Rev. 
Lett. 75, 2466 (1995). 
[12] K. Hara, Y. Sun and T. Mizusaki, Phys. Rev. Lett. 83, 
1922 (1999). 

[13] G.-L. Long and Y. Sun, Phys. Rev. C 63, 021305(R) 
(2001). 

[14] Y. Sun, C.-L. Wu, K. Bhatt, M. Guidry and D.-H. Feng 

Phys. Rev. Lett. 80, 672 (1998). 
[15] Y. Sun, C.-L. Wu, K. Bhatt and M. Guidry, Nucl. Phys. 

A 703, 130 (2002). 
[16] A. Bohr and B.R. Mottelson, Nuclear Structure Vol. I, 

Appendix 1A (Benjanin, Inc., London, 1969). 
[17] Y. Sun and M. Guidry, Phys. Rev. C 52, R2844 (1995). 
[18] Y. Sun, J.-y. Zhang and M. Guidry, Phys. Rev. Lett. 78, 

2321 (1997). 

[19] Y. Sun, J.-y. Zhang, M. Guidry and C.-L. Wu, Phys. Rev. 

Lett. 83, 686 (1999). 
[20] C.-T. Lee, Y. Sun, J.-y. Zhang, M. Guidry and C.-L. Wu, 

Phys. Rev. C 65, 041301 (R) (2002). 
[21] Y. Sun, C.-L. Wu, D.H. Feng, J.L. Egido and M. Guidry, 



Phys. Rev. C 53, 2227 (1996). 
[22] J. A. Sheikh, Y. Sun and P.M. Walker, Phys. Rev. C 57, 
R26 (1998). 

[23] The abnormal-parity (intruder) level in each major shell 
is also taken into account in the FDSM, but treated sep- 
arately. In the context of this paper, we shall ignore it 
for simplicity, since the collective excitations are mainly 
determined by the valence nucleons in the normal-parity 
levels. The detail discussion about the role played by the 
intruder level can be found in the review article of the 
FDSM |I3]. 

[24] For convenience, the definition of P rfi here is twice of that 
in Ref. O- 

[25] A. P. Zuker, J. Retamosa, A. Poves, and E. Caurier, Phys. 
Rev. C 52, R1741 (1995). 

[26] G. Martinez-Pinedo, A.P. Zuker, A. Poves, and E. Cau- 
rier, Phys. Rev. C 55, 187 (1997). 

[27] M. W. Guidry and C. -L. Wu, Int. J. Mod. Phys. E 2, 17 
(1993). 

[28] C.-L. Wu, M.W. Guidry, D.H. Feng and J.-Q. Chen, 

Fizika 22, 123 (1990). 
[29] J.-y. Zhang, MA. Caprio, N.V. Zamfir and R.F. Casten, 

Phys. Rev. C 60, 061304(R) (1999). 
[30] F. Iachello, Phys. Rev. Lett. 87, 052502 (2001). 
[31] Y.S. Chen and Z.C. Gao, Phys. Rev. C 63, 014314 (2001). 
[32] M. Dufour and A.P. Zuker, Phys. Rev. C 54, 1641 (1996). 
[33] T. Bengtsson and I. Ragnarsson, Nucl. Phys. A 436, 14 

(1985). 

[34] J.-y. Zhang, Y. Sun, M. Guidry, L.L. Riedinger and GA. 

Lalazissis, Phys. Rev. C 58, R2663 (1998). 
[35] Y. Sun, J.-y. Zhang, M. Guidry, J. Meng and S. Im, Phys. 

Rev. C 62, 021601(R) (2001). 
[36] P-G Reinhard, Rep. Prog. Phys. 52, 439 (1989). 
[37] P. Ring, Prog. Part. Nucl. Phys. 37, 193 (1996). 
[38] X. L. Han, C.-L. Wu, D.H. Feng and M. W. Guidry, Phys. 

Rev. C 45, 1127 (1992). 
[39] C.-L. Wu, M. Guidry and D.H. Feng, Phys. Lett. B 387, 

449 (1996). 

[40] K. Hara and S. Iwasaki, Nucl. Phys. A 332, 61 (1979); 

348, 200 (1980). 
[41] Y. Sun and K. Hara, Comput. Phys. Comm. 104, 245 

(1997). 



